Universal fluctuations in radial growth models belonging to the 
KPZ universality class 
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PACS 81 . 15 . Aa - Theory and models of film growth 

PACS 05. 40. -a - Fluctuation phenomena, random processes, noise, and Brownian motion 
PACS 89 . 75 . Da - Systems obeying scaling laws 

Abstract. - We investigate the radius distributions (RD) of surfaces obtained with large-scale 
simulations of radial clusters that belong to the KPZ universality class. For all investigated 
models, the RDs are given by the Tracy- Widom distribution of the Gaussian unitary ensemble, in 
agreement with the conjecture of the KPZ universality class for curved surfaces. The quantitative 
agreement was also confirmed by two-point correlation functions asymptotically given by the 
covariance of the Airy2 process. Our simulation results fill the last lacking gap of the conjecture 
that had been recently verified analytically and experimentally. 



Growth phenomena remain a topic of great interest in 
nonequilibrium Statistical Physics, mainly because of the 
self-similarity and universality emerging from dynamical 
local processes of different systems. In this context, one of 
the most important examples is the Kardar-Parisi-Zhang 
(KPZ) universality class introduced by equation [l]: 



dh 



= i/V^/iH- A|V/ip -1-77, 



(1) 



where 77 is a Gaussian noise. This universality class was 
observed in a large number of models [2j[3] and a few ex- 
perimental systems |4||6j. 

Former works on surface dynamics were mainly con- 
cerned with the scaling properties of the surface fluctu- 
ations by means of scaling exponents |2||3|[7|. However, 
there is a number of other universal quantities that are also 
suitable for determining universality in surfaces |6]. Ex- 
amples include the stationary distributions of the global 
interface width and the extremal height (8j|9], and the 
height distributions during the transient regime [lO 11 
that precedes saturation of the interface width. 

The distributions during the transient (growth) regime 
were computed exactly for some models in the KPZ class, 
in 1 -|- 1 dimension 



10 



13| . Among the most relevant cases, 
the height distribution (HD) of the single step model [2], 
determined analytically by Johansson 10 , has the Tracy- 



Widom distribution of the largest eigenvalue of the Gaus- 
sian unitary ensemble (GUE) of the random matrix the- 
as a limit solution 
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Prahofer and Spohn 
also obtained analytically the scaling form of the KPZ 
universality class and a Tracy- Widom distribution for the 
HD of the polynuclear growth model (PNG). Further- 
more, they showed that the HD for the growth from flat 
substrates is given by the Gaussian orthogonal ensemble 
(GOE), while growth from a single seed (radial growth) 
leads to a GUE distribution for the radii fTTlflSl. Re- 



cently, Sasamoto and Spohn 16 found a solution of the 
one-dimensional KPZ equation with an initial condition 
that induces the growth of curved surfaces and the GUE 
distribution was confirmed as the limit for the radius fluc- 
tuations. Subsequently, a numerical evaluation of this so- 
lution corroborated these analytical results 17 . Prahofer 
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and Spohn ^18j also obtained an analytical solution for the 
limiting process describing the surface fluctuations in the 
PNG model as the so-called Airy2 process. 

Experimentally, the GOE and GUE distributions were 
obtained in a few experiments exhibiting KPZ scaling. An 
evidence of the KPZ exponents was found in the slow com- 
bustion of paper sheets |4j . The burning fronts evolve from 
a flat initial condition and the obtained HD has a reason- 
able agreement with GOE. A recent experiment on electro- 
convection of turbulent liquid crystal films [5j allowed to 
investigate an isotropic radial growth with high accuracy. 
It was shown that this system belongs to KPZ universal- 
ity class exhibiting radius distribution (RD) in excellent 
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agreement with GUE, including the cumulants from sec- 
ond to fourth order. However, the mean is shifted and 
tends to the GUE value as a power law The same 



result was found by Sasamoto and Spohn 16 in a solution 



of the KPZ equation in 1 + 1 dimensions with an initial 
edge condition indicating a universal behaviour of the ex- 



ponent 1/3. However, Ferrari and Frings 13 have recently 



shown that the scaling law featuring this approach to GUE 
mean is not universal. Indeed, they have shown that the 
mean shift decays as t~^^^ for the totally asymmetric sim- 
ple exclusion process (ASEP) whereas no correction (up 
order ©(t"^/^)) is found for the weakly ASEP, both mod- 
els belonging to the KPZ class. The liquid crystal film 
setup was also used to induce a front growth from a flat 
surface, and a good agreement with GOE was obtained 
for the HDs g. 

The theoretical and experimental evidences above men- 
tioned strongly suggest that the GUE distribution is an 
universal feature of one-dimensional growth with radial 
symmetry in the KPZ class. However, this conjecture 
is based on a limited number of models with exact re- 
sults 
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and a single experimental work 5|. 
Numerical confirmation of the GUE distributions are, up 
to this moment, missing. In order to fill this gap, we inves- 
tigate the RD of large radial clusters (larger than 3 x 10^ 
particles) generated with different versions of the Eden 
model [19 -21 . We show that the RDs exhibit very good 
agreement with GUE distribution. We also observed that 
the cumulants of the distribution converge to the GUE val- 
ues, as previously observed in other systems (5l|16|. The 



two-point correlation function also converges to the Airy2 



process, as predicted by the conjecture [18 
The Eden model 



19 



consists in adding new particles 
in the empty neighbourhood of a growing cluster. If 
the growth starts with a single particle, the model yields 
asymptotically spherical clusters with a self-afhne surface 
exhibiting the scaling exponents of the KPZ universal- 
ity class 21 . We simulated off-lattice clusters in two- 
dimensions with the usual algorithm [21 22 : a particle 
in the active (growing) zon^ is selected at random and 
a new particle is added in a random position chosen in 
the empty neighbourhood of the selected particle. The 
procedure is repeated while the cluster does not reach N 
particles. With suitable optimizations 
able to grow clusters with up to 3 x 10'' 
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we were 
particles. Since 
we randomly pick up a particle from a constantly updated 
list containing Ng surface sites, the time step is simply 
A< = l/Ng. A total number of up to 10'^ off-lattice clus- 
ters were grown in order to perform statistical averages. 

We also have simulated Eden models on a square lattice 
using an algorithm proposed in Ref. 20 that removes the 
lattice anisotropy effects. The method consists in accept- 
ing a given growth step with probability pj — {rij/AY ^ 
where n„ — 1,2,3,4 is the number of occupied nearest- 
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Fig. 1: (Color online) Radius distributions of on- and ofF-lattice 
Eden models rescaled to a null mean and a unitary variance. 
The mean radius of the aggregates are approximately 2500 for 
the off-lattice model and 3.2 x 10* for lattice models. The 
solid line is the rescaled GUE distribution. In this plot, R* = 
{R-{R})/aR. 



neighbours (NN) of a selected growth site j and is an 
adjustable parameter. This method allows to generate 
isotropic clusters containing more than 4 x 10^ particles. 
We investigated two algorithms for lattice simulations. In 
the version Eden A, one site is randomly selected among 
all Ng growth sites (empty NNs of the cluster) and then 
occupied. In the version Eden B, one of the Nj, sites in 
the cluster border and one of its empty NN are randomly 
chosen and the empty one is occupied. The values of the 
parameter v that produce isotropic clusters are va = 1-72 
and 1/^ = 1 for Eden A and B, respectively [20j. At each 
attempt, the time is increased by l/Ng and 1/iVf, for Eden 
A and B, respectively, independently of the growth suc- 
cess. The numbers of clusters used for statistics were up 
to lO'' for Eden A and up to 10^ for Eden B. Notice that 
statistical fluctuations are much stronger in Eden A due 
to the higher amount of overhangs in the surface. 

We start the analysis with a comparison among the RD 
of the Eden models and the GUE distribution, both suit- 
ably rescaled and shifted to have a null mean and a unitary 
variance. We assume a scaling form 



PiR) ^ 



(2) 



^The active zone consists of tliose particles having sufficient 
empty space in its neigfibourfiood to add at least a new particle. 



where cr|j is the variance of the RD and G{x) is a nor- 
malized scaling function. This scaling form reduces finite 
size corrections, what has improved data collapses in other 
analyses, as interface width distributions [s], for example. 
In Fig. [T] we compare the rescaled RDs for the three Eden 
models with the rescaled GUE distribution. An excellent 
collapse of all curves upon a single curve G{x) was ob- 
tained. Similar results hold for different growth times. 
This data collapse confirms that, a part of corrections to 
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Fig. 2: (Color online) Determination of the parameters related to radius evolution given by Eq. (|3|. (a) Average radius growth 
rate d{R)/dt against t~^^^ for off- (main plot) and on-lattice (insets) Eden models. Dashed lines are extrapolations to f — >■ cxa 
used to determine A. (b) Rescaled interface width for a scale e. The horizontal lines are estimates of the parameter A. Cluster 
average radii are 32000 for on-lattice and 2700 for off-lattice models, (c) The equivalent analysis of (b) for the correlation 
function. 



scaling in the cumulants described below, the RDs of the 
all investigated models agree with the GUE distribution, 
as conjectured by Prahofer and Spohn 
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In radial growth belonging to the KPZ universality 
class, the radii are stochastic variables evolving in time 
as 
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R{t)o,Xt+{A'Xt/2f\GUE 



(3) 



where A and A are two non-universal (model dependent) 
parameters. The random variable xgue is distributed 
according to the GUE Tracy- Widom distribution |14j . 
Therefore, A is the asymptotic radial growth rate obtained 
from A ~ d{R)/dt + ayt'"^'^ in the limit i -J> cx) [23]. In 
Fig. [2f^a) , we show the average radius growth rate against 
a power of time. The extrapolated asymptotic values are 
A « 1.1843(2) for off-lattice, A ~ 0.2639(2) for Eden A, 
and A ~ 0.4807(2) for Eden B, where the uncertainties 
obtained in the regressions are shown in parenthesis. 

The parameter A was estimated in two independent 
ways. We can use the local squared surface roughness, 
in a window of size e, defined as 



w''{e,t)^{[R{x,t)Y),-{R{x,t))l 



(4) 



where (. . .)^ denotes an average within several windows in 
the interface. Alternatively, we can estimate the A value 
using the height-height correlation function 



c(e,i) = {[R{x + e,t)~R{x,t)Y), 



(5) 



For long times, theoretical arguments predict that ~ 
Ae/Q and c c:^ Ae 23 . Curves for 6w^/e as functions of e 



are shown in Fig. [2|^b). Well-defined plateaus are observed 
for both on-lattice models, except for short scales, when 
a small deviation is observed. Since off-lattice simulations 
are much smaller, the plateau is not so evident as in the 



on-lattice case, but the data also tend to a constant value. 
For sake of comparison, we measured the local roughness 
exponent, defined as w{e) ~ e", and found a « 0.43 for our 
largest off-lattice simulations. This value is considerably 
smaller than the expected exponent of the KPZ class a = 
1/2, confirming the presence of strong finite time effects 
in the exponents. Roughness exponent for the on-lattice 
models are very close to the value a = 1/2. The estimates 
of parameter A are represented by dashed lines in Figs. 
[2jb) and (c). The estimates using local interfaces width 
are: Am 1.46 for off-lattice. Aw 1.32 for Eden A, and 
A « 0.71 for Eden B models. The correlation function 
estimates are slightly larger: A w 1.55 for off-lattice and 
A models, and A ss 0.84 for Eden B. 
In agreement to Eq. ([3]), the quantity 

_ R- Xt 

~ (v42At/2)i/3 W 

is a random variable given by a GUE distribution. The 
RDs shown in Fig. |3] were obtained with the values 
A — 1.1842 and A = 1.45 for off-lattice simulations, 
A = 0.263887 and A = 1.43 for Eden A, and A = 0.4806 
and A = 0.805 for Eden B. As predicted by radial KPZ 
conjecture, a very well collapse for different models (and 
different times) is observed. The agreement with the GUE 
distribution is noticeable, except by a shift in q, that van- 
ishes as < — 00. The shift is more evident for Eden B as 
can also be seen in Fig.|4] Our results confirm the limiting 
scaling form conjectured by Phahoffer and Spohn 
as previously observed in theoretical 
experimental [5] systems. 

In order to quantify the agreement between RDs and 
GUE distributions, we investigate the nth order cumulants 
of the probability distribution P{q) denoted by k^j. The 
differences between the cumulants obtained for off-lattice 
simulations and the GUE values are shown, as function of 
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Fig. 3: (Color online) Probability distribution of g = [R ~ 
\t)/{A^\t/2Y''^ for our largest-time simulations; Off-lattice 
Eden model with cluster mean radius f = 2700; Eden A and 
B models with mean radius f = 32000. The solid line is the 
GUE distribution. 

time, in Fig. [ijja). As expected, all cumulants converge to 
the GUE values. As observed experimentally by Takeuchi 
and Sano [s] , the first moment decreases towards the GUE 
value while the higher order cumulants increases towards 
GUE values. The same happens for Eden A, as can be 
seen m Fig. ^h). Differently, the mean for Eden B in- 
creases (more slowly than the others) towards GUE and 
the higher order cumulants converges more quickly to the 
theoretical values. This negative amplitude of the differ- 
ence between simulation and GUE mean was also observed 
in a solution of the KPZ equation 



16 . Indeed, the law 



describing the convergence of the mean have recently at- 
tracted great interest [5 13 16 . Our off-lattice and Eden 
B simulations are very well described by the law pre- 
viously reported [5 16 while the simulation of Eden A are 
only consistent with this approach since a long power law 
regime was not observed. 

A further evidence of the agreement between Eden 
growth and radial KPZ conjecture is yielded by the 
two-point correlation function given by C2(e,i) = 
{R{x + €,t)R{x,t)) — (-R)^, that, in agreement with 
the conjecture [l8], scales at long times as C2{e,t) ~ 
{A^Xt/2)^/'^g2iu) with u = {Ae/2)/{A^Xt/2p\ where 
g2{u) is the covariance of the Airy2 process [24]. In Fig. 
[5] we show the correlation function C2 for different mod- 
els. A very good agreement between the scaled C2 for all 
Eden models and 172 function is obtained, showing that 
the models are well described by the Airy2 process. It is 
worth mentioning that scaled C2 approaches 32 for long 
times and that this approach is slower for on-lattice mod- 
els. 

In summary, the KPZ universality class in radial growth 
has been subject of recent analytical jl^ and experimen- 
tal [5] investigations that agreed with the conjecture pro- 
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Fig. 4: (Color online) Differences between cumulants obtained 
for Eden models and GUE distribution as functions of time 
for (a) off-lattice simulations, (b) Eden A, and (c) Eden B. 
Cumulants up to fourth order are shown in the main plots. 
Insets show the mean against time in a double logarithmic 
scale. Dashed lines have slope -1/3. 
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15 , where interface 



posed by Prahofer and Spohn 
fluctuations in systems belonging to the KPZ universality 
class are described by well-known universal distributions. 
However, a computational verification in growth models 
was lacking until the present work. In the present Letter, 
we have investigated the radius distributions in clusters 
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Airy2 process. We also thank the former discussions with 
Herbert Spohn that motivated the beginning of this work. 
SCF thanks the kind hospitality at the Departament de 
Fisica i Enginyeria Nuclear/UPC. 



Fig. 5; (Color online) Scaled two-point correlation function for 
the Eden models. The average radius of the clusters are 2700 
for off-lattice and 32000 for A and B models, respectively. The 
solid curve is the covariance of the Airy2 process. In this plot, 
C2 = (A2aV2)-2/3c2 and u = {Ae/2){A^ . 



obtained with Eden growth models 20 21 starting from 
a single particle. 

The radius distributions obtained for all models exhibit 
an excellent agreement with the scaling ansatz given by 
Eq. (|3|, that associates the radius fluctuations with the 
Tracy- Widom 14 distribution of the Gaussian Unitary 
Ensemble. The cumulants of order n > 1 associated to 
RD converge to the corresponding GUE cumulants for 
all investigated models. A finite time correction of or- 
der t~^^^ in the first moment was clearly observed for off- 
lattice and Eden B simulations, in agreement with other 
systems [5 16 . Finally, a correlation function in accor- 
dance with the so-called Airy2 process yields a further 
strong evidence that the radius fluctuations in all inves- 
tigated growth models are in agreement with the KPZ 
conjecture. 

As a flnal remark, notice that the the small exponent 
a « 0.43 obtained for the largest off-lattice simulations 
shows that scaling exponents undergo strong flnite time 
effects. Therefore, the RD analysis may be more reli- 
able than scaling exponents to determine the universality 
class of a system. In particular, the scaling form given by 
Eq. ([2]) is simpler than the analysis with Eq. ([s]), since flt 
procedures are not required in the flrst approach and it 
does not have flnite time corrections. 

In conclusion, the Prahofer-Spohn conjecture is now 
fully verified in the three general branches of Statistical 
Physics: experimental, theoretical, and computational. 



This work was partially supported by the Brazilian 
agencies CNPq, FAPEMIG, and CAPES. Authors thank 
F. Bornemann by kindly providing the covariance of the 



REFERENCES 



[1] 
[2] 

[3; 

[4 
[5 

[6; 

[7 

[s; 
[9; 

[fo; 

[11 

[12 
[13 
[14 
[15 
[16 
[17 

[is; 

[19 



[20; 

[21 
[22 
[23; 
[24; 



Kardar M., Parisi G. and Zhang Y.-C, Phys. Rev. 
Lett. , 56 (1986) 889. 

Barabasi A.-L. and Stanley H. E., Fractal Concepts m 
Surface Growth (Cambridge University Press, Cambridge, 
England) f995. 

Meakin p.. Fractals, Scaling and Growth far from Equi- 
libnum (Cambridge University Press, Cambridge, Eng- 
land) 1998. 

L. MiETTiNEN, M. Myllys, J. Merikoski and J. Tim- 
ONEN, Eur. Phys. J. B , 46 (2005) 55. 
Takeuchi K. a. and Sano M., Phys. Rev. Lett. , 104 
(20f0) 23060f. 

Takeuchi K. A., Sano M., Sasamoto T. and Spohn 
H., Set. Rep. , 1 (20ff) 34. 
Krug J., Advances in Physics , 46 (f997) 139. 
Oliveira T. J. and Aarao Reis F. D. A., Phys. Rev. 
E , 76 (2007) 06f60f . 

Oliveira T. J. and Aarao Reis F. D. A., Phys. Rev. 
E , 77 (2008) 04f 605. 

Johansson K., Commun. Math. Phys , 209 (2000) 437. 
Prahofer M. and Spohn H., Phys. Rev. Lett. , 84 
(2000) 4882. 

Kriecherbauer T. and Krug J., J. Phys. A: Math. 
Theor. , 43 (20f0) 40300f . 

Ferrari P. and Frings R. e., J. Stat. Phys. , (20f f) f . 
http : //dx . doi . org/10 . 1007/sl0955-011-0318-4 
Tracy C. and Widom H., Commun. Math. Phys. , 159 
(1994) f5f. 

Prahofer M. and Spohn H., Physica A , 279 (2000) 
342. 

Sasamoto T. and Spohn H., Phys. Rev. Lett. , 104 
(20f0) 230602. 

Prolhac S. and Spohn H., Phys. Rev. E , 84 (20ff) 
0ffff9. 

Prahofer M. and Spohn H., J. Stat. Phys. , 108 (2002) 
f07f. 

Eden M., A two-dimensional growth process in proc. of 
Proceedings of Fourth Berkeley Symposium on Mathemat- 
ics, Statistics, and Probability, edited by Neyman J., 
Vol. 4 (University of California Press, Berkeley,California) 
f96f pp. 223-239. 

Paiva L. R. and Ferreira S. C, J. Phys. A: Math. 
Theor. , 40 (2007) F43. 

Ferreira S. C. and Alves S. G., J. Stat. Mech.: Theor. 
Exp. , 2006 (2006) Pff007. 

Alves S. G., Ferreira S. C. and Martins M. L., Braz. 
J. Phys. , 38 (2008) 8f. 

Krug J., Meakin P. and Halpin-Healy T., Phys. Rev. 
A , 45 (f992) 638. 

Bornemann F., Ferrari P. and Prhofer M., J. Stat. 
Phys. , 133 (2008) 405. 



p-5 



